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We have investigated electronic transport in graphene nanoribbon devices with additional bar-shaped exten¬ 
sions (“wings”) at each side of the device. We find that the Coulomb-blockade dominated transport found in 
conventional ribbons is strongly modified by the presence of the extensions. States localized far away from the 
central ribbon contribute significantly to transport. We discuss these findings within the picture of multiple 
coupled quantum dots. Finally, we compare the experimental results with tight-binding simulations which 
reproduce the experiment both qualitatively and quantitatively. 
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Graphene^ - a monolayer of carbon atoms - pos¬ 
sesses a variety of novel electronic properties including 
a linear energy dispersion relation.^ Nanodevices made 
from graphene were recognized as interesting and poten¬ 
tially useful building blocks for applications and quan¬ 
tum circuits.Narrow graphene stripes called nanorib¬ 
bons®’® are the basic building blocks of more complicated 
graphene nanoelectronic devices. If they have perfect 
edges, they are expected to either exhibit a well defined 
band gap or conducting edge states depending on the 
edge orientation.^ 

Experimentally, most graphene nanoribbons show a 
region of strongly suppressed conductance close to the 
charge neutrality point.At sufficiently low temper¬ 
ature, Coulomb blockade is typically observed in this 
regime of suppressed conductance suggesting that local¬ 
ized charge carriers dominate the electronic transport 
properties.Various mechanisms were suggested to be 
responsible for the observed localization of charge carri- 
ers,5.ii-i6 work has shown that for sufficiently 

clean devices, edge disorder is mainly responsible for 
the experimentally observed suppressed conductance^®’^^ 
and that the sites of localized charge can be close to 
the ribbon edges.In a recent experiment, evidence 
was found, that certain electronic states are mainly local¬ 
ized along the disordered graphene edge and can extend 
along the device edge into the wider graphene leads. 

In this paper we investigate a device geometry which 
we call “nanoribbons with wings”, as the device resem¬ 
bles the outline of a dragonfly. The investigated devices 
consist of a typical nanoribbon design where additional 
graphene stripes of similar width are attached perpendic¬ 
ularly at each nanoribbon edge (see Fig. Id). This device 
geometry was designed to study the interplay between 
transport along the edges and through the bulk of the 
ribbon devices, as the total edge length is significantly in¬ 
creased while the device area connecting source and drain 
contacts is maintained. We show that the additional 
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wings significantly alter the electronic transport prop¬ 
erties of the ribbons. Transport in these winged devices 
is governed by charge localization leading to Coulomb 
blockade. Notably, localized states sitting far out in 
the wings and therefore far away from the direct current 
path contribute to transport. The typical area on which 
charges are localized is estimated using measurements 
in magnetic field. The findings are finally compared to 
tight-binding simulations. 

The two investigated device geometries are shown in 
Figs. Ia,d. The difference between the nanoribbon con¬ 
trol devices (CTRL) shown in Fig. la and the nanorib¬ 
bons with wings (WING) in Fig. Id are the extra 
graphene bars added at each side to the latter ones. The 
two investigated CTRL graphene nanoribbons are about 
220 nm long and 50 nm wide. The central ribbon of the 
three investigated WING devices is of similar length and 
between 50 — 70 nm wide. The added wings are about 
50 nm wide and 550 nm long. All devices were character¬ 
ized in at least two different cooldowns and at different 
applied gate voltages. While details vary, the general 
findings presented in this paper are valid for all investi¬ 
gated devices and cooldowns. 

The mechanically exfoliated^ single layer graphene 
flakes^^’^® were patterned by reactive ion etching.®’^® The 
graphene nanodevices are situated on top of a 285 nm 
thick Si 02 substrate with a global silicon back gate 
(BG). For each device, DG voltages were applied be¬ 
tween source/drain, at each side gate (SG) and at the 
BG. Gurrent between source and drain was recorded as 
a function of the applied voltages. No leakage currents 
could be experimentally detected from any gate to any 
other part of the device. Measurements were performed 
at temperatures of T = 1.6 K and at T = 4.2 K. 

Both GTRL and WING devices were initially charac¬ 
terized by recording the current as a function of the ap¬ 
plied BG voltage at a small constant bias. In agreement 
with previous experiments, a region of suppressed con¬ 
ductance is found.®’®“^® This region is typically at posi¬ 
tive BG voltages due to unintentional chemical doping of 
the devices. The measurements discussed in this paper 
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Figure 1 . (a) Scanning force microscopy (SFM) image of 

one of the control (CTRL) devices. In the brown parts, 
the graphene (yellow) was etched away. The different con¬ 
tacts are marked by black circles. The additional uncon¬ 
tacted graphene stripes above and below the ribbon were pat¬ 
terned in order to keep the device geometry similar to the 
WING devices, (b) Current flowing between S and D in one 
CTRL device as a function of the applied side gate voltages 
(VsD = 1 mV, T = 4.2K). (c) Current flowing through one 
CTRL device as a function of the applied back gate voltage 
and source-drain bias {T — 4.2K, VsG 3 = VsG 4 = OV). (d) 
SFM image of one of the ribbons with wings (WING): ad¬ 
ditional graphene stripes were patterned at the side of the 
main graphene nanoribbon, (e) Current flowing through the 
WING device as a function of the applied side gate (SG) volt¬ 
ages (T = 1.6 K, VsD = ImV). (f) Current flowing through 
the WING device as a function of the applied BG and bias 
voltages (T = 1.6 K, VsGi = Vsg2 = OV). For the measure¬ 
ments in (b,e), the same ranges in SG voltage were applied. 
For (c,f), the same range in BG voltage was applied, whereas 
a different range in bias voltage was used. In (c,f), the current 
is plotted on a logarithmic scale. 


are all recorded in this regime of suppressed conductance. 

Transport through the CTRL devices is governed by 
Coulomb blockade (CB) in agreement with previous stud¬ 
ies.When measuring the current as a function of the 
applied BG and bias voltage at constant applied SG volt¬ 
ages, CB diamonds are observed as shown in Fig. Ic. 
Further information about the position of the localized 
charges resulting in CB is obtained by following the po¬ 
sition of a CB peak while changing the applied gate volt¬ 
ages.The current flowing through one CTRL device 
as a function of applied SG voltages is shown in Fig. lb: 
diagonal lines with a slope of minus one indicate that the 
localized states couple equally well to both SGs. This is 
not surprising as the two symmetric SGs are located far 
away from the ribbon and therefore cannot discriminate 
potentially existing localized states at different positions 
inside the ribbon. 

Transport in the WING devices is also governed by 
GB as shown in Fig. If. However, the diamonds are 
about a factor of 10 narrower in BG and a factor of 4 
smaller in bias direction. Also, the evolution of GB peaks 
in SG voltages is clearly different as shown in Fig. le: 
many closely spaced and oscillating lines are found run¬ 
ning nearly vertically or horizontally through the mea¬ 
surement. Additionally, a modulation of the current am¬ 
plitude is found when varying the SG voltages. 

The nearly vertical/horizontal resonances observed in 
the WING devices show that areas much closer to one 
SG than to the other have an influence on transport. 
This would, for example, not be the case in a diffusive 
and incoherent metallic system. In one of the WING de¬ 
vices, the outermost 200 nm of the wings were removed by 
RIE after the device was characterized. This new device 
with shortened wings exhibited qualitatively and quanti¬ 
tatively similar transport properties to the unshortened 
version. The only notable exception is that the slopes 
of the nearly vertical/horizontal lines got flatter/steeper, 
indicating that the outermost 200 nm of the wings did 
contribute to the transport characteristics. 

Further information is obtained by measuring the cur¬ 
rent as a function of the applied BG and SG voltages as 
shown in Figs. 2a,b: both plots look very similar at a first 
glance, despite the fact that in one case SGI (a) and in 
the other case SG2 (b) was used. This indicates that the 
part of the device being equally strongly coupled to both 
SGs, i.e. the central ribbon part, determines the over¬ 
all conductance. The finer lines with flatter slopes (in¬ 
dicating stronger coupling to the SGs) differ for the two 
measurements in Figs. 2a,b, indicating that the transport 
details are also influenced by the wings. 

Next, the spatial extent of localized states is investi¬ 
gated. Fig. 2c shows transport data of a WING device 
obtained in a perpendicular magnetic field. Such plots 
were also recorded as a function of the SG voltages and 
show similar behavior. Up to magnetic fields of about 
1 T, the resonances remain mostly unchanged. In the re¬ 
gion of 1 — 2 T, the GB resonances start to move, broaden 
and increase in amplitude. At higher fields, starting 
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Figure 2. (a,b) Current flowing through one WING device as 
a function of one SG and the BG voltage. The other SG is 
kept at ground potential (T = 1.6 K, Vsd = ImV). While the 
overall behavior is similar, the fine features are distinctly dif¬ 
ferent. As a guide to the eye, two features with different slopes 
are marked with black lines, (c) Current flowing through one 
WING device as a function of the applied BG voltage and the 
perpendicular magnetic field (T = 1.6 K, Vsgi = Vsg2 = OV, 
Vsd = ImV). The current stays unchanged up to fields of 
1 —2 T. At fields of about 3 T, the current increases by roughly 
one order of magnitude (note logarithmic colorscale). 


from 2 — 3 T, the current suddenly increases by more 
than an order of magnitude. Assuming that the addi¬ 
tion of one flux quantum to the area of a localized state 
will significantly change transport^^^^®, the following ar¬ 
eas can be extracted: A(1T) = h/{eB) « (65 nm)^ and 
A(2T) « (45 nm)^. The strong increase in current above 
2 — 3 T can be attributed to reduced backscattering in 
electronic transport due to the chirality induced by the 
magnetic field.Alternatively, one could try to estimate 
the area of a localized state from the width of a CB dia¬ 
mond. This approach is however not useful here as it is 
experimentally challenging to determine which resonance 
belongs to which localized state. 

In summary, the experimental findings suggest a pic¬ 
ture of tunneling coupled states of varying spatial ex¬ 
tent whose locations are spread throughout the whole 
structure including the remote ends of the wings. These 
states interact via Coulomb interaction^^’^^, resulting in 
a huge multi-quantum-dot molecule. Due to the overlap 
of the wave-functions, current flow through the structure 
can also be influenced by states localized far away from 
the central ribbon. The enhanced number of Coulomb- 
blockade diamonds in Fig. If is the manifestation of an 
enhaced number of states per energy window and rib¬ 
bon length due to the device area being increased by the 
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Figure 3. Current flowing through a 1.6 pm long CTRL-type 
device as a function of the applied BG and bias voltage (T = 
4.2K, no side gates). The size of the CB diamonds is about 
an order of magnitude larger than for the WING devices. 


wings. In such a situation, the distinction between cur¬ 
rent flowing along the edge or within the bulk becomes 
meaningless. This contrasts with a CTRL-type ribbon 
of increased length shown in Fig. 3, where the area is 
increased but the number of states per energy window 
and ribbon length remains the same. As a result, the 
diamonds overlap strongly due to the serial arrangement 
of the localized states. 

In order to confirm this general picture and to obtain 
more insights into the spatial distribution of wave func¬ 
tions, we have performed tight-binding simulations using 
third-order nearest neighbor tight binding.For an anal¬ 
ysis of localized eigenstates, we consider a central ribbon 
element with one set of wings of experimental dimen¬ 
sions (see Fig. 4). We connect the top and bottom of 
the center region using periodic boundary conditions to 
obtain a simulation cell of manageable size. We assume 
no bulk disorder and a randomly jagged edge with an 
edge roughness of the order of 2 nm (with 5 nm correla¬ 
tion length). We identify four different types of wave- 
functions as shown in Figs. 4a-d: (a) states mostly delo¬ 
calized over the whole structure; (b) states mostly located 
in the central ribbon; (c) “square” states localized in the 
wings, coupling both edges; and (d) states localized pre¬ 
dominantly along one edge of the wings. The delocalized 
states (a) are mostly found far away from the Dirac point. 
When approaching the Dirac point, fewer states couple 
directly through the central ribbon (b). Instead, states 
localize either at a single edge (d), or between the two 
edges of the wing (c). 

To extract length scales from the simulations, we cal¬ 
culated scattering states inside a wing geometry. We find 
exponential (Anderson) localization in direction x along 
the wing 

= J y)f dyKi Ax exp(-x/Lioc). (1) 

The localization length Lioc(e) as a function of energy is 
extracted by averaging over 80 scattering states in a nar¬ 
row energy window (e ± 5meV). Tioc(e) decreases close 
to the Dirac point (see red curve in Fig. 4e). A local¬ 
ization length of the order of 50 nm (the ribbon width) 
is consistent with experiment. States located in either 
flap contribute to Coulomb blockade peaks, where the 
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Figure 4. Envelopes of states simulated for the geometry of a 
central ribbon part with one set of attached wings, (a) State 
mostly delocalized, (b) State coupling directly through the 
middle ribbon element, (c) “Square” state coupling both rib¬ 
bon edges together, (d) A state mostly localized along one 
ribbon edge, (e) Averaged localization length Lioc of scatter¬ 
ing states along the wing as function of energy [see Eq. (1)], 
for three different magnetic field strengths. Each square rep¬ 
resents an average over 80 scattering states. 


overlap of their exponential tails with the center region 
of the dot is responsible for the tunneling coupling into 
and out of the localized state. Indeed, we numerically 
find identical localization behavior when directly prob¬ 
ing the decay of eigenstates away from their maximum 
(not shown): “square” states decay exponentially away 
from their center region (of width ss 50 nm) due to edge 
roughness. 

At finite magnetic fields of T ss IT, the localiza¬ 
tion length slightly increases as shown in Fig. 4e (pur¬ 
ple curve). As soon as twice the magnetic length As ~ 
25nm/v'5/T is smaller than the ribbon width, Lioc 
abruptly increases, in good agreement with the experi¬ 
mental findings in Fig. 2c. 

In summary, we have presented electronic transport 
experiments of a special graphene nanoribbon geometry 
with added graphene wings at each side of the ribbon. 
While transport close to the charge neutrality point is 
still governed by Coulomb blockade, the details change 
strongly. We show that even the outer parts of the wings 
contribute to transport. From transport measurements 
in perpendicular magnetic field, we extract a typical lo¬ 
calization area being smaller than the wing size. Previ¬ 
ous experiments have observed states localized along the 
edge.^^ Such states are again found by the tight-binding 
calculations performed for the ribbon with wings de¬ 
vices. The comparison of the ribbons with wings with the 
long control ribbons shows however, that a picture where 


transport is only governed by charges localized along the 
edge is insufficient to explain the measurements. Instead 
we propose a picture where the different localized states 
are tunneling coupled. In such a picture it is meaning¬ 
less to distinguish transport along the edge or over the 
bulk. This picture is compatible with both experiment 
and tight binding simulations. Notably, the simulations 
manage to reproduce both the typical localization length 
extracted from the experiment as well as the experimen¬ 
tal behaviour in a magnetic field. 
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